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ABSTRACT 

We study the evolution of galaxy structure since z ~ 1 to the present. From a Great Observatories 
Origins Deel Survey South (GOODS-S) multi-band catalog we define (blue) luminosity- and mass- 
weighted samples, limited by Mb < —20 and M* > 10 10 M©, comprising 1122 and 987 galaxies, 
respectively. We extract early- type (ET; E/SO/Sa) and late- type (LT; Sb-Irr) subsamples by their 
position in the concentration-asymmetry plane, in which galaxies exhibit a clear bimodality. We find 
that the ET fraction, /et, rises with cosmic time, with a corresponding decrease in the LT fraction, 
/ltj m both luminosity- and mass-selected samples. However, the evolution of the comoving number 
density is very different: the decrease in the total number density of Mb < — 20 galaxies since z = 1 
is due to the decrease in the LT population, which accounts for ~ 75% of the total star formation rate 
in the range under study, while the increase in the total number density of M* > 10 10 M galaxies 
in the same redshift range is due to the evolution of ETs. This suggests that we need a structural 
transformation between LT galaxies that form stars actively and ET galaxies in which the stellar 
mass is located. Comparing the observed evolution with the gas-rich major merger rate in GOODS-S, 
we infer that only —20% of the new ET galaxies with M* > 10 10 M appeared since z ~ 1 can be 
explained by this kind of mergers, suggesting that minor mergers and secular processes may be the 
driving mechanisms of the structural evolution of intermediate-mass (M* - 4 x 10 10 M©) galaxies 
since z — 1. 

Subject headings: galaxiesievolution - galaxies interactions - galaxies structure 



1. INTRODUCTION 

Local galaxies present two main populations in the 
color- magnitude diagram: the red sequence, formed pri- 
marily by old, spheroidal quiescent galaxies, and the blue 
cloud, formed p rima rily by spiral star- forming galaxies 
(jStrateva et al-ll2QQlb iBaldrv et all 120041 ). It is now well 
establi shed that such b imodality is p resent at higher red- 
shifts dBell et all 12004 up to z ~ 1; Cass ata et all 120081 : 
lllbert et al.l I2010L up to z ~ 2; Kriek et al. 2008L at 



These results are not immediately expected from 
the popular hierarchical A-CDM models, in which 
the more massive dark matter halos are the final 
stage of successive mergers of smaller halos. How- 
ever, whether downsizing poses an important problem 
for A-CDM depends on the correct understanding of 
the baryonic phy sics, for which models are making re- 
cent progress (seelBower et al. 2006; D e Lucia fc Blaizotl 



z — 2.3), and appears to be strongly linked to mass: 
more massive galaxies were the first to finish forming 
their stars and populating the red sequenc e. This mass 
depen dence has been dubbed downsizing ([Cowie et al.l 
1996): massive galaxies having experienced most of their 
star formation at early times and being passive by z — 1, 
and many among the less massive gala xies experience ex- 
tended star-fo r mation histories (see iBundv et aLl 120061 : 
iScarlata et all l2QQ7al [P erez- Gonzalez et al.l |2008, and 
references therein). 
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120071 : iStewart et al.l 120091 : iHopkins et alJl2QQ9dL and ref- 
erences therein), as well as on proper accounting of the 
dependency of the halo merger history with environment. 
Because of that, the role of galaxy mergers in the buildup 
of the red sequence, and their impact on the evolution of 
galaxy properties, i.e., color, mass, or structure, remains 
an important open question. 

In addition, the redshift evolution of the mass function 
suggests that the red sequence grows because star for- 
mation is quenched in the blue cloud (|Bell et al.l 120071 : 
Ruhland et al. 2009). Because -80% of red sequence 
galaxies ar e morphological early types (ETs) at z < 1 
(E/SQ/Sa. IStrateva et all 120011 : iLotz et al.ll2008bf ). and 
because the star formation is loca ted in spiral galaxies 
(jBell et al.ll2005Ujogee et al.ll2QQ9h . the blue-to-red tran- 
sition may be accompanied by a late- to ET transforma- 
tion. 

In this paper we study the role of gas-rich major merg- 
ers, selected by morphological criteria, in the structural 
evolution of galaxies since z — 1. In a previous paper, 
lLopez-Sanjuan et al.l (|2009ai hereafter L09), we study 
the major merger rate evolution in GOODS-S, finding 
that only - 10% of z = galaxies with M* > 10 10 M 
have undergone a gas-rich major merger since z ~ 1. 
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Our goal is, first to quantify the structural evolution, 
and, second, to compare it against the merger rate. 

The paper is organized as follows: in Section [2] we sum- 
marize the GOODS-S data set used, and in Section [3] we 
develop the methodology to determine the fractions of 
ET and late- type (LT) galaxies versus redshift. Then, 
in Section |4] we summarize the obtained early- and LT 
fractions and comoving number densities, and their evo- 
lution with z, while in Section [5] we study the role of 
gas-rich major mergers in the observed ET evolution. 
Finally, we present our conclusions in Section [6] We 
use H = 70 km s -1 Mpc~\ Vt M 0.3, and Q A 0.7 
throughout. All magnitudes are Vega unless noted oth- 
erwise. 

2. GOODS-S DATA SET 

2.1. Galaxy samples 

We work with the galaxy catalog from the Great 
Observatories Origins Deep Survey S outh (GOODS- 
S) 7 field by the Spitzer Legacy Team (jGiavalisco et al.l 
2004). We used the Version 1.0 catalogs 8 and re- 
duced mosaics in the (£435), F606W (V 606 ), 
F775W (z 775 ), and FSbOLP (z 850 ) Hubble Space Tele- 
scope/Advanced Camera for Surveys {HST/ ACS) bands. 
These catalogs were cross-correlated using a 1.5" search 
radius with the GOODS-S Infrared Array Camera 
(IRAC) selected sampl e in the Rainbow co s molo gical 
data base 9 published in | Perez-Gonzalez et al.l (|2008l see 
also iPerez-Gonzalez et al.ll2005l and G. Barro et al., in 
preparation), which provided us spectral energy distri- 
butions (SEDs) in the UV-to-MIR range, well-calibrated 
and with reliable photometric redshifts, stellar masses, 
star formation rates, and rest-frame absolute magni- 
tudes. 

We refer the reader to the mentioned papers for a 
more detailed description of the data included in the 
SEDs and the analysis procedure. Here, we summarize 
briefly the main characteristics of the data set. The Rain- 
bow database contains consistent aperture photometry 
in several UV, optic al, NIR, and MIR bands w ith the 
method described in Per ez-Gonzalez et al.l (|2008f ). UV- 
to-MIR SEDs were built for 4927 IRAC sources in the 
GOODS-S region down to a 75% completeness magni- 
tude [3.6] =23.5 mag (AB). These SEDs were fitted to 
stellar population and dust emission models to obtain 
estimates of the photometric redshift (z p hot)> the stellar 
mass (M*), and the rest-frame B-band absolute magni- 
tude (M B ). 

Rest-frame absolute £?-band magnitudes were esti- 
mated for each source by convolving the templates fit- 
ting the SED with the transmission curve of a typical 
Bessel-5 filter, taking into account the redshift of each 
source. This procedure provided us accurate interpo- 
lated 5-band magnitudes including a robustly estimated 
if-correction. Stellar masses were estimated using the 
exponential sta r formation PEGASE01 models with a 
Salpeter (1955) initial mass function (IMF) and various 
ages, metallicities and dust contents included. The typ- 
ical uncertainties in the stellar masses are a factor of 

7 http:/ /www. stsci.edu/science/goods/ 

8 http://archive.stsci.edu/prepds/goods/ 

9 http: / /guaix. fls.ucm.es/~pgperez /Proyectos / 
ucmcsdatabase.en.html 



~2, typical of most stellar population st udies (see, e.g. 
iPapovich et"aIll2QQa IFontana et al.ll2006f ). 

The median accuracy of the photometric redshifts at 
z < 1.5 is |^ S pec — %>hot|/(l + Z snec) = 0-4, with a fraction 
<5% of catastrophic outliers (jPerez-Gonzalez et al.l l2008. 
Figure B2). In the present paper we use a Zphot = o~s z (1 + 
%>hot) as %>hot uncertainty, where <js z is the standard 
deviation in the distribution of the variable S z = (%>hot — 
2spec)/(l + 2phot)> that is well des cribed by a Gaussian 
with mean fi$ z ~ and as z (see Lopez-Saniuan et al. 
I2009U for details). We take a 6z = 0.043 for z < 0.9 
sources and o~s z = 0.05 for z > 0.9 sources. 

From the Rainbow catalog described above we defined 
two samples in the range 0.1 < z < 1.3. One 
sampl e is selected in lum inosity, Mb < —20 (~ M B at 
z ~ 0; iFaber et al.l 120071 ). which comprises 1122 sources. 
The value of Mr is 1 mag brighter at z <y 1 than lo- 
cally (jGabasch et al.ll2004l : lllbert et al.ll2005l : IFaber et al.l 
l2QQ7h . evolution that reflects the descent in the star for- 
mation rate density of t he universe since z ~ 1 (e.g. 
iHopkins fc Beaco m 2006). Hence, if we assume a con- 
stant Mb cut, we select different areas of the luminosity 
function at each redshift, biasing our results. However, 
L09 show that asymmetry as reliable morphological in- 
dicator in GOODS-S is valid only for M B < -20 galax- 
ies, so we decided to use this constant cut, instead of 
an evolving one, to ensure good statistics. This selection 
condition our results, so is important take it into account 
when we interpret the evolution of the comoving number 
density of galaxies with redshift (Section l47Tj) . The sec- 
ond sampl e is selected in mass, M+ > 10 1Q M (~ 0.1M* 
at z ~ 0; IPerez-Gonzalez et al.l 12008). which comprises 
987 galaxies. In this case M * evolves little, if any, since 
z ~ 1 dPozzetti et al.l 120071 : IPerez-Gonzalez efaLl 120081 : 
iMarchesini et al.l l20Q9h . so we decide to use a constant 
selection in mass that ensures 75% completeness for pas- 
sively evolving galaxies in the range under study (see 
L09, for details). 

2.2. Morphological Indices 

We use conc entration (C; lAbraham et al.lll994f ) and 
asymmetry (A: lAbraham et al.lll996[ ) indices to perform 
our structural study. Concentration is defined as 



C = 5 x Ioq 



^80 

r20 



(1) 



where r2o and rgo are the circular radii which contain 
20% and 80% of the total galaxy flux, respectively. Con- 
centration correlates with several properties of galaxies, 
as bu lge-to-total ratio ([Strateva et al. 2001; Conselice 
2003), absol ute ff-band mag nitude (|Conselicd l2QQ3h . or 
stellar mass ([Cons elice 2006a). For details about concen- 
tration measurements see fBershad v et all |2000). 
The asymmetry index is defined as 



A = 



£14)-/: 
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where Io and Bo are the original galaxy and background 
images, Jiso and Biso are the original galaxy and back- 
ground images rotated 180°, and the summation spans 
all the pixels of the images. The background image is a 
sky source- free section of 50 x 50 pixels. This index gives 
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us information over the source distortions and, we can 
use it to identify recent m erger s ystems which are highly 
distorted (e.g.. iConselicd l2QQ3t iDe Propris et all l2QQ7t 
iBridge et all 120071 : iLopez-Saniuan et al.ll2009bD ~All de- 
tails about the measurement of the asymmetry index in 
the HST/ ACS images of GOODS-S sources are given in 
L09, with similar techniques having been applied for the 
measurement of the C index. To avoid statistical mor- 
phological ^-corrections and to deal with the loss of in- 
formation with redshift (i.e., spatial resolution degrada- 
tion and cosmo logical dimming) , we determined C and A 
indices in £?-band rest-frame galaxy images, which were 
artificially redshifted to a unique, representative redshift, 
Zd = 1. This provides us an homogeneous data set to 
perform structural studies in GOODS-S. 

2.3. Major merger rates 

We use the major merger rates from L09 thought this 
paper. In that work, we selected as major merger rem- 
nants those galaxies with high values of the asymme- 
try index (A > 0.3). We performed that study for 
M* > 10 10 Mq galaxies up to z ~ 1 in the same mass- 
selected sample that we use in the present paper (Sec- 
tion 12. ip , and we took into account the effect of observa- 
tional errors in z and A, which overestimate the merger 
fraction due to the spillover of normal sources to the high- 
asymmetry regime, by maximum likelihood (M L) tech- 
niques developed in lLopez-Sanjuan et al.l ((2008). We ob- 
tain lower merger fractions (/™ ph ^$5%) tha n other de- 
terminations also based on m orphology (e.g., lLotz et al.l 
2008b; Conselice et al. 2009), but are in good agreement 
with those from 1 Lopez-Sanju an et al.l (|2009b[ ) in Groth 
Strip apply ing the same metho dology as that of ours, and 
those from Uogee et al.l j2009) in Galaxy Evolution from 
Morphology and SEDs 10 (GEMS) by eye-ball inspection 
of the sources. To obtain t he merger rate (jft m ) fr o m the 
merger fraction, we assume IPerez-Gonzalez et al.l ([2008) 
mass functions and a typical timescale of Ta ~ 0.5 Gyr 
(|Conselicell2006bL 120091 : lLotz et al.ll2008ai 12009( 1. 

The methodology in L09 is only sensitive to gas-rich 
major mergers (i.e., at least one of the merging galaxies 
is a LT, gas-rich galaxy): asymmetry as merger indicator 
is calibra ted in local univ erse with ultraluminous infrared 
galaxies ([Con selice 2003), while Ta is determin ed by TV- 
body s imulations of gas- rich major mergers ([Conselicel 
2006 rj lLotz et al . 2008a). In a forthcoming paper, we 
show that nearly all the high- asymmetry galaxies in our 
sample are also star- forming galaxies, supporting the no- 
tion that we are only sensitive to gas-rich mergers (C. 
Lopez-Sanjuan et al, in preparation). 

3. THE CONCENTRATION-ASYMMETRY BIMODALITY 

Concentration and asymmetry indices are useful 
to segregate galaxies by their structure, i.e., bulge- 
dominated galaxies (E/SO/Sa, ET galaxies in the fol- 
lowing), and disk-dominated galax ies and irregulars (Sb- 
lrr, LT galaxies in t he following; Bersha dv et al.l [2 000: 
Lauger et al. 2005a; Menanteau et al. 2006; Yagi et al. 
20061 : iHuertas-Companv et al.l f2008: N eichel et al.ll2008D . 
To segregate galaxies by t heir morphology we need an in- 
dex as t he cumpliness (S; IConselicd 120031 ), or the bumpi- 
ness (B; Blakeslee et al. 2006), which give us information 

10 http: / /www. mpia.de / GEMS /gems. htm 
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Fig. 1. — Distribution in the concentration-asymmetry plane of 
the 757 galaxies in the catalog with Mb < — 20 and 0.35 < z < 1.1. 
The black lines mark the different structural areas: ET galaxies 
(E/SO/Sa), LT galaxies (Sb-lrr), and gas-rich major mergers. The 
contours show number of galaxies (see the right legend). 

about the distribu tion of the star formation in the galaxy 
(|van der W ei 2008). On the other hand, high values of 
the asymmetry index select gas-rich major merger rem- 
nants (Section 12. 3p . 

In Figured] we show the distribution of the 757 galax- 
ies in the catalog with Mb < —20 and 0.35 < z < 1.1, 
whic h present two d i fferen t po pulations, in agreemen t 
with iZamojski etaD (|2QQ7h and iConselice et~atl (|2QQ9D . 
One population is associated with ET galaxies, with its 
maximum at (C,A) = (3.25,0.09), and the other with 
LT galaxies, with its max i mum a t (C,A) = (2.52,0.17), 
as shown by lllbert et al.l (|2006l ). We determined the 
positions of the two peaks in (C, A) space as follows: 
first we fitted the histogram in concentration space with 
two Gaussians, which gave us the concentration values 
of the two peaks as well as the minimum between them 
(C = 2.9). Then, we fitted the asymmetry histogram 
of C < 2.9 and C > 2.9 galaxies with a simple Gaus- 
sian to obtain the asymmetry maxima. In Figure [T] 
we also show the three areas that we use to segregate 
galaxies: gas-rich major mergers (A > 0.3, see L09), ET 
galaxies (A < 0.3 n A < 0.3C - 0.75), and LT galaxies 
(A < 0.3 fl A > 0.3C - 0.75). The limit A = 0.3C - 0.75 
was chosen to optimize the ET/LT separation (see the 
following section for details). 

3.1. Early- and Late-Type Fraction Determination 

The first goal of this paper is to study the fraction 
of ET (/et) and LT galaxies (/lt) as a function of 
redshift. To obtain reliable ET and LT fracti ons we 
used ML techniques. lLopez-Sanjuan et al.l (j2008[ ) devel- 
oped a two-dimensional ML method to determine merger 
fractions from the asymmetry in galaxy images; the 
method is far superior to classical source counting, as, 
with source counting, measurement errors cause source 
spillover to neighboring bins. The ML method has been 
successfully used to measure the gas-rich merger frac- 
tion in Groth strip (Lo pez-Sanjuan et al . 2009b) and 
in GOODS-S (L09), and we apply it here for the de- 
termination of /et and /l t- We refer the reader to 
lLopez-Sanjuan et al.l (|2008f ) for all the details about the 
ML method and its assumptions. The main difficulty in 
the application of the ML method in the present study 
is that it works in two dimensions, whereas our problem 
has three dimensions: C, A, and the redshift z. 

We could work in (C, z) space, and use as selection 
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TABLE 1 

Early- type fraction in GOODS-S at 0.6 < z < 0.85 



Sample Selection 


™LSS 


w/ LSS a 


w/o LSS b 


LSS (z = 0.735) c 


M B < -20 
M* > 10 10 M 


72 
94 


Q07+0.029 

u.oy / _o.027 
u - ouy -0.027 


n oo 9 + 0.033 
UlOOZ -0 .030 
n 4QQ+0-035 


6 o 4 +0.045 
u.uy / _ 047 



a Early-type fraction in the sample with LSS sources. 
b Early-type fraction in the sample without LSS sources. 
c Early- type fraction in the LSS. 



limit the condition C = 2.9, which marks the separation 
between ET and LT populations in concentration space 
(see the previous section). However, this selects galaxies 
with A > 0.3, i.e., gas-rich major mergers, as ET and LT 
galaxies. Indeed, the separation between ETs and LTs is 
oblique in Figure [T] and is well described by A = 0.3C — 
0.75 (drawn in the figure). If we define a second variable 
C A = A — 0.3C + 0.75, LT galaxies and major mergers 
have CA > 0, while ET galaxies have CA < 0. That is, 
the fraction of ET galaxies is /et = fcA<o, while the 
fraction of LT galaxies is /lt = fcA>o — Zm ph > where 
f™ ph is the morphological merger fraction reported by 
L09. The error in the new variable is Gq A — a\ + 
0.09a^. 

How are the results affected by the selection limit? We 
determined /et with three different limits: C = 2.9, 
A = 0.3C-0.75, and A = |C-2.30. We find that all the 
values of /et are consistent within la with our preferred 
condition, A — 0.3C — 0.75, in the three redshift ranges 
under study, namely, z\ = [0.35, 0.6), z 2 = [0.6,0.85), 
and zs — [0.85, 1.1). We select these particular ranges to 
resemble those in L09. Because the particular election of 
the selection limit does not bias our results, we use the 
limit A = 0.3C - 0.75 in the following. 

3.2. Large- Scale Structure Effect 

It is well known that the more prominent large-scale 
structure (LSS) in the GOODS-S field i s located at red- 
shift z = 0.735 (Raviku mar et al.ll2007D . The next two 
more important ones are located at z — 0.66 and z — 1.1. 
The former is an over density in redshift space, but not in 
the sky plane, while the latter is a cluster, but comprises 
an order of magnitude fe wer sources than t he z = 0.735 
structure (145 versus 

i2,Ed 

ami et al.H2005r >. Because of 
this, we concentrate on the LSS at z = 0.735 and ignore 
other structures in GOODS-S. In order to check the effect 
of this LSS in our derived structural fractions, we recalcu- 
lated the ET fraction in the range Z2 = [0.6, 0.85) by ex- 
cluding the sources within Sv < 1 500 km s -1 (Sz ~ 0.01) 
of z = 0.735 (jRawat et al.l 120081 : L09). In Table [Q we 
summarize the number of sources in the LSS for each 
sample (riLss)? the ET fractions in the LSS and the ET 
fractions in the field, both in the samples with and with- 
out LSS sources. The value of /et in the samples with- 
out LSS sources is ~0.07 lower than the determination 
in the whole samples. More importantly, both values 
are incompatible at la. This fact exp lains the high ET 
fraction found by lLauger et al.l (j2005b) at that redshift 
in their GOODS-S field study. Because of this, in the 
following we use the structural fractions obtained from 
the samples without LSS sources. On the other hand, 
/et S ~ ^S/et^ that is > the ET fraction is higher in 
dense environments. In Figure [2j we show the position 
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Fig. 2. — Position in the sky plane of the M* > 10 10 M galax- 
ies in the more prominent LSS in GOODS-S field (z — 0.735). 
Red dots are ET galaxies, blue squares are LT galaxies, and green 
triangles are gas-rich major mergers. 



TABLE 2 

Structural fractions of 
M b < -20 GALAXIES 



z 


/et 


/lt 


0.475 
0.725 
0.975 


n o 71 +0.057 
U -^' i -0.048 
n oqo+0.033 
UlOOZ -0 .030 
n 99P ;+0.027 

u.zzo_ 024 


603+ - 051 

U.OU«3_q 059 
6 o 7 +0.032 
u - uo ' -0.037 
n 7 o 9 +0.026 
u - ' oz -0.031 



in the sky plane of ET (red dots), LT (blue squares), and 
major mergers (green triangles) for M* > 10 10 M galax- 
ies in the LSS: ET galaxies are concentrated in more pop- 
ulated regions, while LT galaxies are located in the out- 
skirts of the stru cture, as expe cted by the morphology- 
density relation (|Dresslerlll980[ ). In addition, the three 
gas-rich major mergers in the LSS are also located in 
the outskirts, a natural consequence of the external lo- 
cation of their progen itors, i.e. LT galaxies (see also 
iHeiderman et all I2009D . 

4. THE STRUCTURAL EVOLUTION OF GOODS-S 
GALAXIES 

We summarize in Table [2] the structural fractions for 
Mb < —20 galaxies in three redshift ranges: z\ — 
[0.35, 0.6), z 2 = [0.6, 0.85), and z 3 = [0.85, 1.1). The frac- 
tion of ET galaxies increases with cosmic time, while LT 
galaxies fraction decreases. Note that errors are formal 
and do not t ake into account c osmic variance, denoted 
a v ; following iSomerville et al.l (j2004f ) , we expect a v ~ 
20%. In the top panel of Figure [3l we show our results 
(squares) with those from the literat ure: diamonds are 
the ET and LT fractions provide bv lLotz et al.l (|2008bl ) 
for Mb < —18.8 — 1.32 galaxies in All- Wavelength 
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Structural fractions of 



n 1 1 1 r 

/et" E-Sa /lt ■ Sb-lrr 




M* > 10 10 M G 



0.0 0.2 0.4 0.6 0.8 1.0 

Z 

Fig. 3. — ET (red) and LT galaxy fraction (blue) vs. red- 
shift. Top: st ructural fractions for Mg < —20 galaxies : this 
work (squares), Lotz et all (I2008U diamo nds), lllbert et al.l (2006, 
invert ed triangles), Scarlata et al. ( 2007b, triangles), IDriver et al.l 
(2006, stars), and Conselice (2006a, pentagons). Bottom: struc- 
tural fract ions for M* > 10 1 Mq galaxies: this work (circles), 
and Mandelbaum et al. (2006, triangles). In both panels, the red- 
shift errors show the redshift range covered for each data, while 
red (blue) line is the linear least-squares fit to the /et (/lt) data. 



Extended Groth Strip International Survey (AEGIS 11 ; 
I Davis et al.l [2007), which mimic our select i on crit eria 
at z ~ 0.9. Pentagons are from Conselice (2006a]) for 
M B < -20 galaxies in th e RC3 12 (Third Reference C ata- 
logue of Bright Galaxies: Ide Vaucouleurs et al.l li991): we 
take ET = E/S0 and early-disks, while LT = late-disks 
and irregulars (his Table 1). Triangles and inverted trian- 
gles are obtained integrating the morphological luminos - 
ity functions (M LF) provided bv | Scarlata et all (j2007bl ) 
at z = 0.7, and lllbert et all (|2006f ) at 0.05 < z < 1.2, 
respectively. In both cases, we integrated the MLF for 
galaxies brighter than M B = -20. lllbert et al.l (|2006f ) 
perform their study in the GOODS-S field and segregate 
galaxies in ET/LT by their positio ns in the C — - A plane , 
a similar methodology to ours. Scar lata et al.l (j2007bl ) 
use Zurich Estimator of Structural Types (ZEST) in Cos- 
mological Evolution Survey (COSMOS 13 ; Scoville et"a!1 
|2007[ ) to classify ~10,000 galaxies in several morpholog- 
ical types: E/S0 (Tl), disks (T2, divided in four sub- 
types by their Sersic index) and irregulars (T3). We 
take ET = Tl + T2.0, while LT = T2.1 + T2.2 + T2.3 
+ T3 -/™ ph (0.7), where /™ ph (0.7) is the morphologi- 
cal merger fraction at z = 0.7 from L09. This selects 
disks wit h the Sersic index n > 2.5 as bulge- dom i nated 
galaxies (Mandel baum et all 120061 : iTruiillo et all 120071 : 
iDahlen et al.ll2007l ). Finally, stars are obtained integrat- 

11 http://aegis.ucolick.org/ 

12 http://heasarc.nasa.gov/W3Browse/all/rc3.html 

13 http:/ /cosmos. astro. caltech.edu/index. html 





> 10 10 Mq 


GALAXIES 


z 


/et 


/lt 


0.475 
0.725 
0.975 


547+0.055 
439+ - 035 

u ^ oy -0 .033 

349+ - 034 
u.o^y_ 031 


445+ - 051 

U.^O_ Q 056 
con+0.034 

u.ooy_ 037 

614 +0-033 
u,D14 -0 .038 



ing the n > 2 (E Ts) and n < 2 (LTs ) luminosity func- 
tions provided by IDriver et al.l (j2006[ ) from Millennium 
Galaxy Catalog (MGC 14 : lEiske et alJl2003l ). They found 
that their z ~ 0.1 galaxies present a bimodality in n, 
with a minimum at n ~ 2. In the three previous cases 
we determine the uncertainty in structural fractions by 
varying lcr, the value of Mg in the integration. The top 
panel of Figure [3] shows that, when we use a similar lu- 
minosity and structural selection, all /et and /lt are in 
good agreement. The linear least-squares fit to the data 
yields 



/etW = (0.56 ±0.02) - 
/lt(*) = (0-43 ±0.02) 



- (0.35±0.02)z, 

- (0.30 ±0.03)2. 



(3) 
(4) 



With these fits, we infer that ET galaxies are the domi- 
nant Mb < — 20 population by number since z ~ 0.2. 

In Table [3] we summarize the structural fractions for 
M* > 10 10 M Q galaxies. The fraction of ET galaxies 
increases with cosmic time, as in the previous case, but 
/et in the mass-selected sample is higher than in the 
luminosity-selected sample in all the redshift intervals 
under study. In the bottom panel o f Figure |3] we show 
our r esults (circles) with those from Man delbaum et all 
(|2006l . triangles). They provide / ET (n > 2.5) and / LT 
(n < 2.5) for -33000 M* > 10 1 Q M galaxies from Sloan 
Digital Sky Survey (SDSS 15 ; Adelman- McCarthy etHI 
2006). The linear least -squares fit to the data yields 



/et(*) = (0.60 ±0.04)- 
fur(z) = (0.40 ±0.04) 



(0.24±0.06)z, 
- (0.19±0.06)z. 



(5) 
(6) 



With these fits, we infer that ET galaxies are the dom- 
inant M > 10 10 Mq population by the number since 
z — 0.5, a higher redshift than in the luminosity-selected 
sample. Although in the range of study a linear function 
is a good approximation to the observed evolution, we 
expect that a power-law function may provide a better 
para meterization when h igher redshift data are available 
(e.g. iFontana et aDl2009D . 

4.1. Comoving Number Density Evolution 

To better understand the structural evolution of galax- 
ies since z = 1, in this section we study the comoving 
number density of ET (pet) and LT galaxies (pur) as a 
function of redshift: 

p ET (z, M) = f ET (z, M)p tot (z, M), (7) 
p LT (z, M) = f hT ( z , M)p tot (z, M), (8) 

where M = Mb [M*] is the selection criteria, and 
ptot(z : M) is the comoving number density at redshift 

14 http:/ /www. eso.org/~jliske/mgc/ 

15 http://www.sdss.org/ 
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Fig. 4. — ET (red) and LT comoving number density (blue) 
vs. redshift. Top: structural number den s ity for Mb < —20 
galaxies: this w ork (squares), Lotz et al. (2008b, diamonds), 
Ilbert et al. (2006, inverted triangles), Scarlata et al. (2007b, trian- 
gles), Driver et al. (2006, stars), and Conselice (2006a, pentagons). 
Bottom: stru ctural fractions for M* > 1 10 Mq galaxies: this work 
(circles), and Mandelbaum et al. (2006, triangles). In both panels 
the redshift error bars show the redshift range covered for each 
data, while red (blue) line is the linear least-squares fit to the pet 
(pur) data. 



TABLE 4 

Comoving number densities of M b < 
GOODS-S 



-20 GALAXIES IN 



z 


PET 


PLT 


ptot(M s ) 




(10- 3 Mpc- 3 ) 


(10- 3 Mpc" 3 ) 


(10- 3 Mpc" 3 ) 


0.475 


1.02 ±0.15 


1.73 ±0.15 


2.77 


0.725 


1.16 ±0.10 


2.26 ±0.12 


3.50 


0.975 


0.96 ±0.11 


3.16 ±0.12 


4.28 



z of galaxies more luminous [massive! than Mb JM±L To 
obtain ptot(z, Mb), we assume the lFaber etHI (|2QQ7D lu- 
minosity functi on parameters, while to obtain ptot(^? M*) 
we assume the iPerez-Gonzalez et al.l ([20081 ) mass func- 
tion parameters. 

We summarize in Table 2] the comoving number den- 
sity for Mb < — 20 galaxies. The number density of LT 
galaxies decreases with cosmic time, while ET number 
density is roughly constant. In the top panel of Figure HJ 
we show the number density data of this work and those 
from literature, obtaining by applying Equation ([7]) and 
dEJ) to the fractions in Figure 03 The linear least-squares 
fit to all data yields 



lo gio(PET 
logioOur) 



) = (" 



3.01 ±0.03) - 
-3.14 ±0.02) 



(0.03±0.04)z, (9) 
-(0.67 ±0.02)*. (10) 



The number density of LT galaxies at z = is 
the value at z = 1, while pet ~ constant. 



20% 



TABLE 5 

Comoving number densities of M* > 10 10 M 

GALAXIES IN GOODS-S 



z 


PET 


PLT 


ptot(M*) 




(10- 3 Mpc" 3 ) 


(10- 3 Mpc" 3 ) 


(10- 3 Mpc" 3 ) 


0.475 


3.41 ±0.33 


2.82 ±0.36 


6.29 


0.725 


2.12 ±0.16 


2.60 ±0.17 


4.83 


0.975 


1.28 ±0.12 


2.29 ±0.13 


3.71 



The total number density of Mb < — 20 galaxies de- 
creases by a factor of 3 since z = 1, a reflection of 
the global decline in the star forma tion rate density 

of the universe as ~ [1 ± z) 4 (e.g. iLillv et al.l [1996; 

Hopkins & Beacom 2006; Perez-Gonzalez et al. 2005; 



Tresse et al. 2007r iVillar et all 120081 ). Because both the 



number densities of LT galaxies and that of gas-rich ma- 
jor mergers decrease with cosmic time (e.g. L09), we ask 
which of these two populations drives the decrease in the 
star formation rate density of the universe. 

We define the ratio between the star formation rate in 
a given population (SFR pop ) and the star formation rate 
in major merger systems (SFR m ) as: 



R 



SFR 



pop 



Mpop/j 



pop 



pop 



SFR n 



mph J^ 1 



(ii) 



where f™ ph and / pop are the morphological major merger 
fraction (i.e., gas-rich, Section I2T3]) . and the fraction of 
galaxies of a given population, respectively, /i pop is the 
fraction of galaxies in the population that are forming 
stars actively, e$F is the star formation in an interacting 
system respect to the median in isolated galaxies, T e is 
the time over which the star formation is enhanced in 
an interaction, and T m ^ is the timescale over which one 
major merger is selected as high-asymmetric source by 
L09 methodology. In Equation (fTT|) . we have assumed 
that all the gas-rich major mergers have enhanced star 
formation, and we do not take into account dry mergers, 
that have negligible star formation. With this definition, 
the total star formation rate (SFR to t) in a given redshift 
range is 

SFR tot = (Rut ± Ret + l)SFR m , (12) 

while the fraction of the total star formation located in 
major merger systems is 

/8F ' m = Jfcr + kr + r ^ 

Becau se most of the red galaxies are ETs (jLotz et al.l 
2008b), we assume that fi^T — (i-e., all ET galax- 
ies are passive), and Ret = 0. On the other hand , 
most of the blue galaxies are LTs (jLotz et al.l [2008b). 
so we take jutt = 1. With these assumptions the star 
formation rate is more important in LT galaxies than 
in major merger systems if Rjjt > 1 and vice versa if 
Rlt < 1- To compare the results with previous works, 
we restrict our study to the range 0.4 < z < 0.8, in which 

/-P h = 0.026tnSg and /lt = 0.6101^32- We assume 
that T m . A = 0.475 ± 0.125 Gyr (L09), esp = 1-50 ± 0.25 
(iRobaina etHI I2009L see als o iLin et al.l 120071 : iLi et al.l 
120081: iKnaoen fc J ames 2009) , and T ( = 2.0 ± 0.25 Gyr 
(jdi Matteo et al.ll2007l : ICox et all feOOS). With these val- 
ues we obtain Rut = 3.7^2.5? that is, the bulk of the star 
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formation in the range 0.4 < z < 0.8 is located in LT 
galaxies, and major mergers account for ~20% of the to- 
tal star formati o n rate, in quanti t ative agreement w ith 
Bell et all l200gjWolf et al.l I2005L Uogee et al.l 1200! or 



Sobral et al.1j2 009). 



Following Roba ina~et al.l (|2009f ) , we can also obtain the 
fraction of SFR to t that is triggered directly by mergers: 



.rtri 
/SF,m 



esF 



"/SF,] 



esF - 1 



esF(#LT + Ret + 1) 



(14) 



If we assume that Ret — 0, as previously, then 
fsF m = ^—2%? a l° w value in good agreement with 
IRobaina et all (|2009l ): they infer that f$ m = 8 ± 3% 

in GEMS 16 (|Rix et al.l 12004ft , at the same redshift 
range, and by studying the correlation function of 
> 10 10 Mq blue galaxies. If we repeat this study 
in the range 0.85 < z < 1.1, we obtain Rut = 2.4 (i.e., 
~ 30% of the total star formation is located in major 
mergers at z ~ 1), and /g g m ~ 9%. Our values are also 
in good agreement with iHopkins fc Hernquistl j2009); 
they inferred the merger-induced burst history up to 
z ~ 2 from the luminosity profiles of local E/S0 galaxies, 
finding that ~ 5%-10% of total star formation is 
triggered by mergers, independently of redshift. Finally, 
our results are the same if we use the mass-selected 
sample. 

In Table [5j we summarize the number densities for 
M* > 10 10 Mq galaxies. In this case the behavior is 
the opposite than in the luminosity-selected sample: the 
number densities of ET and LT galaxies increase with 
cosmic time. This increase is more important for ET 
galaxies. In the bottom panel of Figure HJ we show the 
number density dat a of thi s work (circles) and those from 
iMandelbaum et ahl ([2006, triangles). The linear least- 
squares fits to the data are 

logio(pET) = (-2.19 ± 0.04) - (0.68 ± 0.06)*, (15) 
logio(pLT) = (-2.37 ± 0.03) - (0.29 ± 0.05)*. (16) 

The ET population increases its number density by a 
factor of 5 between z = 1 and * = 0, while LT galaxies 
increase by a factor of 2 in the same range. Because 
~40% of M* > 10 10 galaxies in the local uni verse 
are in place at * = 1 ([Perez-Gonzalez et al.l 12008). the 
increase in the number density of these galaxies since 
* = 1 is primarily due to the increase in pet- As in the 
previous section, we expect that a power-law function 
may provide a better parameterization th an the linear 
when higher redshift data are available (e.g. jTavlor et all 
2009). If we fit a power-law function to the pet data, 
Pet (2) oc (1 + z) a , we obtain a = —2.3 ± 0.2, while 
iTavlor et al.l (|2009h find a = -1.70 ± 0.14 for M* > 
10 11 Mq red sequence galaxies. If we assume that most 
red sequence galaxies are ETs (jLotz et al.l 1200 8b). this 
implies that red massive galaxies evolve less in number 
densi ty than less m a ssive ones since * ~ 1, in agreement 
with lFerreras et al.l ([2009D . 

The results of this section suggest that we need a 
structural transformation in the range < * < 1 be- 
tween LT galaxies, which form stars actively, and ET 

16 http://www.mpia-hd.mpg.de/GEMS/gems.htm 



galaxies, in which stellar m ass is located, in agree- 
ment with (Bell et al.l (|2007f ). IVergani et all (|2008f ). or 
iRuhland et al.l (j2009[ ). The question is, can gas-rich ma- 
jor mergers drive this structural transformation? 

5. THE ROLE OF GAS-RICH MAJOR MERGERS IN THE 
STRUCTURAL EVOLUTION OF INTERMEDIATE-MASS 
GALAXIES 

To explore the role of gas-rich major mergers in the 
structural evolution of M* > 10 10 M galaxies we define 
the fraction of new ETs that appears between Z2 and z\ 
due to gas-rich major mergers, /ET,m? as: 



/ET,m(^l,^) 



(21,22) 

p^ e T v (z 1 ,Z 2 ) PET (2i) - PET (22) 



(17) 

where Pet, m (21, ^2) is the number density of new ETs 
due to gas-rich mergers, and Pet^i,^) is the total 
number density of new ETs. We assume th at each gas- 
rich major merger remnant is an ET galaxy (jNaab et al.l 
2006al iRothberg fc JosepHl2006al lbl: IHopkins et al.l l2008. 
2009b|), so pET,m (21,22) = Prem (21,2:2), being 



Prem(2l,Z 2 ) = / ' 5R m (0) (1 + z)^ 1 - 
J z-i 



dz 



H E(zY 



(18) 



where E(z) = y^^a + ^m(1 + z) 3 in a flat universe, 
3?m(0) = (0.3 ± 0.1) x 10- 4 Mpc- 3 Gyr- 1 , and n = 
3.5 ± 0.4 (merger rate parameters are those from L09). 
Using the results of the previous section, we infer that 
/ET,m(0, 1) = 17^7°%. The errors take into account the 
uncertainties in the merger rate and number density pa- 
rameters. This value implies that gas-rich major mergers 
cannot explain the observed structural evolution. 

Interestingly, if we extrapolate the observed tendencies 
up to Z2 = 1.5, we obtain that /et,hi(1, 1-5) ~ 100%. 
This is, gas-rich major mergers could be the dominant 
process in ETs formation at z J> 1. Although extrapo- 
lating our results beyond z = 1 is risky and more data in 
the range 1 < z < 2 are needed to explore the suggested 
picture, the very different behavior in the two redshift 
intervals makes our results qualitatively reliable. 

How are our results modified by spheroid-spheroid 
(dry) mergers? The L09 meth odology is only sensitive to 
gas-rich mergers (Section [22]), so the role of dry mergers 
cannot be measured directly in this work. However, we 
can estimate the effect of dry mergers in our results: if we 
assume that a dry merge r between two ET galaxies leads 
to another ET galaxy ([G onzalez- Garcia & van A lbadal 
2003: iGonzalez-Garcia et all 120061: iNaab et al.l l2QQ6bt 
iBovlan-Kolchin et al.l 120061 : IHopkins et al.l l2009d h then 
each dry merger remnant represents the disappearance 
of one ET galaxy. That is, the value of Pet^i,^) in 
Equation (fT7|) could be higher due to dry mergers. In 
addition, some of our observed gas-rich mergers may 
be spiral-spheroid (mixed) mergers that conserve the 
number of ET galaxies (i.e., only spiral-spiral wet merg- 
ers create new spheroids). Because of this, the value of 
Prem(2i, Z2) in Equation ([17]) could be lower due to mixed 
mergers. In summary, our /ET,m(0, 1) value is at least an 
upper limit to the importance of gas-rich major merg- 
ers in the structural evolution of M* > 10 10 M galaxies. 
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5.1. Comparison with Previous Studies 

In this section, we compar e our result with those in lit- 
erature. Bu ndv et al.l (|2QQ9l ) perform a close pair study 
in GOODS-South and North, and compare the forma- 
tion rate of new spheroids against their measured gas- 
rich (wet + mixed) merger rate. From their results, we 
infer / ET ,m(0.4, 0.9) ~ 16±?g% for log(M*/M ) ~ 10.7 
galaxies 17 that compares well with our value in the same 
redshift range, /ET,m(0.4, 0.9) = 26^0%, f° r the mam 
galaxy in our mass sample, which has log (M*/M ) ~ 
10.6. Although we include Sa galaxies in our definition 
of ET galaxy, the conclusion from both works is simi- 
lar: less than half of the new intermediate-mass ET that 
appeared at 0.4 < z < 0.9 come from gas-rich major 
mergers. 

In their work. lWild et al.l(j2009[ ) studied post-starburst 
(PSB) galaxies in COSMOS and their role in red se- 
quence assembly. The spectra of PSB indicate that the 
formation of O- and early-B-type stars has suddenly 
ce ased in the galaxy, wh ile the simulations performed 
by Uohansson et al.l ([20081 ) find that the PSB phase can 
only be reached by gas-rich major merger remnants. In 
L09, we show that the gas-rich merger rate and the PSB 
rate of M* > 10 10 M galaxies at 0.5 < z < 1.0 are 
similar, supporting PSB galaxies as descendants of our 
gas-rich mergers. In addition, most of the red sequence 
galaxies are ETs (e.g. llLotz et al.ll2008b[ ). so we can com- 
pare the fraction of red sequence mass from PSB galaxies 
at 0.5 < z < 1.0, 38t^%, with / ET ,m(0.5, 1) = 34±^%. 
Despite the different methodologies and assumptions, the 
agreement between both studies is remarkable. 

Other w orks also study the impact of mergers in ETs 
evolution. iBundv et al.l ((2007) compare the evolution of 
the virial mass functions of halos hosting spheroidals with 
the merger rates predicted by cosmological simulations, 
concluding that major mergers are insufficient to ex- 
plain th e observed increase in ET population. lLotz et all 
(2008b) use morphological indices to study the major 
merger rate to z ~ 1.2 in a B— band selected sample 
(Lb > 0AL* B ) and compare it against the evolution of 
E/SO/Sa galaxies. They find that all the ET evolution 
can be explained by major mergers, although the effect 
of observational errors that tend to overestimate sys- 
tematically the merger fra ction by morphological indices 
([Lopez- Sanjuan et al.ll2QQ8h makes their result an upper 
limit. 

Summarizing, our results are in good agreement with 
previous works when a similar mass selection is ap- 
plied. At higher stellar masses the picture is dif- 
ferent: the merger fraction depends o n stellar mass 
(jde Ravel et al.l 120091 : IBundv et al1l2009[ ). having M* > 
10 11 Mq galaxies higher pair fractions than less massive 
ones, and being red pairs more common at these masses 
(Bun dv et al.1120091 ). This suggests dry mergers as an im- 
portant pro cess in the evolution of massive galaxies since 
z ~ 1 (e.g. iBell et "all 120041: lEinet al.l 12008: Il bert et al.l 
20101). The size ([T ruiillo et al. 2007; Buitrago et al.l 



2008; Ivan Dokkum et al.l 120081 : Ivan der Wei et all 
and velocity dispersion evolution (jCenarro fe Tr uiillo 
l2QQ9h of M* > 10 11 M ET galaxies since z ~ 

For our cosmology and assuming a Salpeter (1955) IMF. 



2 also support the importance of mergers, specially 
the impact of minor mergers in the evolution of 
these systems (Bezanson eTaT1 l200l iNaab et al.l [2009; 
iHopkins et al. 2009a). This problem was also analyzed 
bv lEliche- M oral et al.l ([2009), who modeled the evolu- 
tion of luminosity function backwards in time for M* > 
10 11 Mq galaxies, selected according to their colors 
(red/blue/total) and their morphologies. They find that 
the observed luminosity function evolution can be natu- 
rally explained by the observed gas-rich and dry major 
merger rates, and that 50%-60% of today's E/S0 in this 
mass range were formed by major mergers at 0.8 < z < 1, 
with a small number evolution since z = 0-8 (s ee also 
iCristobal-Hornillos et all 120091 : lllbert et all 12010ft . Note 
that the gas-rich major merger fractions assumed by 
lEliche-Moral et al.l (|2009ft are those from L09 for 5-band 
selected galaxies (Mb < — 20), which were obtained in 
a similar way as the merger fractions used through this 
paper. 

Th is makes M* - 10 11 M© (iPerez-Gonzalez et al.l 
2008) a transition mass at z < 1: at higher masses major 
mergers are an important process in the evolution of ET 
galaxies, while at lower masses other mechanisms dom- 
inate the observed evolution (see also iDrorv fc Alvarezl 
l2CMIyajTder Wei et al.1l200l lOescn et al.ll2009ft ~ 

5.2. Which Processes are Responsible of Early-type 

Rise? 

Since gas-rich major mergers are not the dominant pro- 
cess in the late- to early-type transformation at interme- 
diate masses, two possibilities remains: minor mergers 
and secular processes. We measure the structure in the 
5-band rest frame, so the increase in the Sersic index 
that we observe can be either the result of an increase 
in the mass of the spheroidal component of the galaxies, 
or due to the decline of the star formation in the disk 
of late types. In the form er, minor mergers increas e the 
Sersic index of galaxies (jEliche- Moral et a l. 2006) and 
multiple minor mergers can l ead to ETs (jBournaud et al.l 
[2003; iHo^k ins et al.l 12009 oft , while secular evolution, 
e.g. bars and disk instabilities, produces pseudob- 
ulges (iKormendv fc Kennicuttl 120041 : iFisher et all 120091 : 
Combes 2009). In the latter, minor mergers can mod- 
ify the gas distr ibution of th e galaxy, shutting down the 
star formation ([Bekki 1998), whereas gas exhaustion or 
the prevention of star formation due to the stabilization 
of the gas in a bulge-dominat ed galaxy (morphological 
quenching, iMartig et al.l 12009) is an example of a secu- 
lar process. To obtain new clues about ET evolution, 
we study the star formation properties of ET and LT 
galaxies in GOODS-S in a forthcoming paper (C. Lopez- 
Sanjuan et al., in preparation). 

6. CONCLUSIONS 

We have studied the structure of Mb < —20 and 
M* > 10 10 M© galaxies in GOODS-S at z < 1. We use 
the position of galaxies in the concentration-asymmetry 
plane to segregate them in ET galaxies (bulge dominated, 
E/SO/Sa), LT galaxies (disk dominated and irregular, 
Sb-Irr), and gas-rich major mergers (see L09, for details). 
We find that: 

1. The ET fraction increases with cosmic time in both 
luminosity- and mass-selected samples, while the 
LT fractions decrease; 
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2. The number density of LT Mb < — 20 galaxies de- 
creases with cosmic time, while ET number density 
is roughly constant. We infer that star formation is 
located primarily in late spirals and irregulars, in- 
stead that in major merger systems. These systems 
account for a ~ 20%— 30% of the total star forma- 
tion rate in the range 0.4 < z < 1.1, while the star 
formation triggered directly by interactions in this 
range is f^ m - 8%; 

3. The number density of ET M* > 10 10 M© galaxies 
increases by a factor of 5 since z ~ 1. This im- 
plies that we need a structural transformation in 
the range < z < 1 between LT galaxies, which 
form stars actively, and ET galaxies, in which stel- 
lar mass is located; 

4. When we compare the observed structural evolu- 
tion with the gas-rich merger rate in the range 
< z < 1, we obtain that only ~20% of the newly 
formed ET galaxies that appear since z = 1 can 
be gas-rich major merger remnants, whereas in the 
range 1 < z < 1.5 these mergers can explain all 
the inferred structural evolution. This suggests 
minor mergers and secular processes as the prin- 
cipal mechanisms in the rise of intermediate-mass 
(M* - 4 x 10 10 Mq) ET galaxies since z ~ 1, 
while gas-rich major mergers may be the domi- 
nant process at higher redshifts (z>l) and masses 
(M* > 10 11 M ). 

We study the star formation properties of GOODS-S 
galaxies in a forthcoming paper (C. Lopez-Sanjuan et 
al., in preparation) to obtain new clues about ETs rise 
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